Collective oscillation in two-dimensional fluid 
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Large-scale collective oscillation is discovered in the two-dimensional Euler equations. For initial 
conditions far from a base stationary flow, the system does not relax to the base stationary flow, but 
instead shows pairs of coherent vortices moving along the base stream line, which leads to large-scale 
oscillatory fields. The investigation of the vicinity of a bifurcation point suggests that this oscillation 
appears through Ifopf bifurcation. Furthermore, a dynamic self-consistent theory explains that this 
oscillation results from the collective organization of a state of self-oscillation. 
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Two-dimensional (2D), and similarly geophysical, 
flows of fluids with high Reynolds numbers are widely 
observed in nature. Examples include oceanic and at- 
mospheric flows [1]. 2D flows are also matters of exper- 
iment [2]. A characteristic phenomenon in such flows 
is large-scale pattern formation. Jupiter's red spot and 
Kuroshio of the north pacific ocean are well-known ex- 
amples in nature. Although highly turbulent on inter- 
mediate and small scales, stationary patterns of flow ap- 
pear on large scales. This is due to the so-called inverse 
cascade [3], which is an essential difference from three- 
dimensional turbulence. 

A theoretical interest lies in the statistical mechanics 
description of the large-scale patterns [4, 5], which has 
been studied intensively over the last two decades [6]. 
Instead of the infinite resolution description of the bare 
vorticity field, a probabilistic description of the field on 
a coarse-grained space is introduced. That is, the proba- 
bility density function of the vorticity level is defined at 
each point. Then maximizing the corresponding Shannon 
entropy functional yields the large-scale stationary pat- 
tern of flow. With the use of this theory, Jupiter's red 
spot, for example, has been successfully explained [7]. 

Here we recall the general scheme of macroscopic sys- 
tems. Needless to say, equilibrium states are described 
with statistical mechanics as entropy maximum states. 
When the systems are non-equilibrium but close to equi- 
librium, they show simple relaxation to, or fluctuation 
around, the equilibrium states. When the systems are 
far from equilibrium, they often show temporal motions 
on macroscopic scales, such as periodic oscillation, by 
spontaneously breaking the temporal symmetry [8] . 

In this context, we reflect the 2D fiuids. The large- 
scale stationary patterns are regarded as the equilib- 
rium states of statistical mechanics, as mentioned above. 
When the systems are slightly far from the stationary 
flows, they show relaxation to [9, 10], or fluctuation 
around [11, 14], the stationary flows, though their be- 
haviors are not so simple as ordinary thermodynamic 
systems. Then, by thinking parallelly to the above route 
to non-equilibrium, we expect large-scale temporal mo- 
tions also in the 2D fluids. Indeed, it has been reported 



that macroscopic oscillation appears, through Hopf bi- 
furcation, in the 2D Navier-Stokes equations under ex- 
ternal shear force [12]. However, this shear obviously 
breaks the symmetry, which may cause temporal mo- 
tions. Then, as a nontrivial question, it is natural to 
ask if the large-scale non-stationary motions exist with- 
out net external forcing. The present paper reports the 
discovery of such large-scale temporal oscillation in the 
2D Euler equations. 

We consider the 2D Euler equations, 

dtuj + v-Vuj = 0, (1) 

on a doubly periodic domain V — [— L2./2, Lj./2) x 
[—Ly/2,Ly/2). Here v = {vx,Vy) is the velocity fleld, 
and Lo — dxVy — dyVx is the vorticity fleld. We set = 27r 
and Ly = 2ttT, with the aspect ratio F > 1, without loss 
of generality. 

The large-scale motion is realized when the initial con- 
dition is set far from the stationary state obtained from 
the statistical mechanics description. We give initial con- 
ditions in the form. 



(2) 



where l^(r) is a stationary base flow and (5a;(0,r) is an 
initial perturbation. 

As the base flow, we consider the so-called zonal flow, 
or the Kolmogorov flow, of the form. 



27r 27r 
n{r) = -— cos— y 

L,y Ly 



(3) 



It is easily shown that this is a stationary solution of 
the 2D Euler equations. Moreover, this is the entropy 
maximum state in the statistical mechanics description 
of the 2D Euler equations on a doubly periodic domain 
[13, 14]. 

Since we are interested in large-scale behaviors, it is 
natural to apply the initial perturbations on the largest 
scale. We consider the initial perturbations of the form. 



27r 27r 
oujyt = 0, r) = cos —x 



(4) 
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where e represents the distance of the initial condition 
from the stationary state, which is not necessarily small. 

In the numerical simulations, we add a very weak 
hyper-viscous term (— l)^+^i^V^'*aj to the r.h.s. of (1), 
with h = 4 and z/ = 2 • 10~^*, just to stabilize the nu- 
merical scheme. We use the pseudo-spectral method [15], 
which is known as the most precise and robust numer- 
ical algorithm, and is indeed regarded as the standard 
method, for the 2D Euler equations in doubly periodic 
domains. We use the resolution 256 x 256. 

After a short transient (t < 200), the system goes into 
stable states. We find three sorts of states, depending 
on the parameters (e, T). To quantatively distinguish 
these states, we introduce an order parameter that is the 
Fourier component of the smallest wavenumber. 



(5) 



where Ldq{t) = uj(t, r)e Note that the imagi- 

nary part is kept zero for the present initial conditions. 

One of the states is stationary zonal flow (FIG. l-(a)), 
similarly to the base flow on the largest scales. Accord- 
ingly, Z{t) is almost zero, with small fiuctuations (FIG. 2- 
(a)). Another state is stationary dipole fiow (FIG. l-(b)); 
the vorticity field shows a dipole pattern on the largest 
scales. Accordingly, Z{t) is almost kept to a nonzero 
value, with small fluctuations (FIG. 2-(b)). Note that 
these two patterns are relevant to the entropy maximum 
states in the statistical mechanics description [13, 14]. 
These mean that, on the largest scales, the system re- 
laxes (strictly, close) to the entropy maximum state, just 
as an ordinary thermodynamic system relaxes to equilib- 
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FIG. 1: Snapshots of the relaxed vorticity field; (a) zonal 
(r = 1.5, e = 0.1) and (b) dipole (F = 1.1, e = 0.5) flow. 

The other state, and what we are interested in, is 
the non-stationary pairwise motion of coherent vortices 
(FIG. 3- (a)). Two pairs of positive and negative coher- 
ent vortices point-symmetrically move along the stream 
line of the base fiow. The motion of these vortices leads 



to oscillatory field. Indeed, Z{t) shows oscillation with a 
rather large amplitude (FIG. 2-(b)). This is in a strong 
contrast to the above two states in that no relevant state 
exists in the statistical mechanics description. 
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FIG. 2: Typical time series of the order parameter Z{t), for 
(a) (r,e) = (1.5,0.1) (zonal), (b) (1.1,0.5) (dipole), and (c) 
(1.5,0.5) (oscillatory). 
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FIG. 3: (a) A snapshot of the vorticity field in the 2D Euler 
equations (F = 1.5, e = 0.5). and (b) the Poincare section in 
the one-body dynamics (7) {b = 0.042, T = 12), when the 
phase of oscillation is zero. F — 1.5. 

By changing the parameters e and F, we obtain the 
phase diagram, which is shown in FIG. 4. 

The system shows a transition from the stationary 
zonal flow to the oscillatory flow. Then we investigate 
how this transition occurs. FIG. 5 shows the amplitude 
of oscillation of Z{t) against the intensity of perturba- 
tion e. The amplitude shows a continuous bifurcation 
at a critical point. In particular, the amplitude increases 
proportionally to (e — Ec)^^^, where ec is the critical point. 
This result suggests that the bifurcation from the zonal 
to oscillatory flow is a Hopf-type bifurcation. 

Next, we investigate the mechanism how the large-scale 
oscillation is kept stable, by a dynamic self-consistent 
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FIG. 4: The phase diagram. 
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FIG. 5: The bifurcation diagram. The amphtudc of oscilla- 
tion oi Z{t) against e and e — Ec (inset in log- log plot), with the 
curve proportional to (e — Ec)^^^, where we obtain Cc = 0.25. 
r = 1.5. The error bars denote standard deviations of the 
small fluctuations. 



theory [16]. Suppose an ensemble of independent "test 
point vortices" moving along the stream line of a tempo- 
rary oscillating auxiliary field on the largest scale. The 
collective motion of the point vortices results in a dy- 
namic field in the limit of infinite number of the vortices. 
If these two fields, one being the cause and the other the 
effect, are identical, this system is no different from the 
original autonomous 2D Euler equations. 

Suppose, on V, the temporally periodic auxiliary vor- 
ticity field on the largest scale in the form. 



'(t) cos — 



(6) 



where Z^'"^{t) is the order parameter of the auxiliary field 
being sinusoidal with respect to time with the amplitude 
h and period T, = bsm2Trt/T. Each of the test 

point vortex moves along its stream line, following the 
equations of motion, 



X = u^^'^ = +dyilj'''''' = sin ^ 



(7a) 



r Ottt 

y = n^"^ = -5, V""" = -2Z-""(i) sin — (7b) 



where ipaux is the corresponding auxiliary stream func- 
tion, w'^"'^ = -V^t/)*"''. 

We consider this one-body dynamical system with h 
and T the same value as Z{t) numerically obtained from 
the original 2D Euler equations. FIG. 3-(b) shows the 
Poincare plot of the trajectory on 2-Kt/T = (mod 27r), 
which is similar to the standard map [17]. In between 
the region of fixed point {y = 0,Ly/2) and Kolmogorov- 
Arnold-Moser tori {y = ±Ly/4), there are the so-called 
chaotic seas. In addition, within the chaotic seas, we find 
two pairs of large 1-to-l resonance island of tori. 

We compare in FIG. 3 this Poincare section of 
the one-body dynamics (b) with the snapshot of the 
vorticity field of the original 2D Euler equations (a) 
when the phase of oscillation is each zero. We no- 
tably find that the position of the cohcirciiit vortices 
in the 2D Euler equations coincides to that of the 
1:1 resonance island of tori in the one-body dynam- 
ics; {x/L^,y/Ly) = (±0.25, tO.125), (±0.25, ±0.625) at 
27Tt/T = (mod 27r). 

This coincidence self-consistently explains the mecha- 
nism of the large-scale oscillation. Once a considerable 
amount of vorticity is concentrated on the 1:1 resonance 
islands of tori, it leads to large-scale oscillatory field. 
Thus generated oscillatory field, in turn, holds the vor- 
ticity inside the islands, and prevents it from diffusing 
away. Thereby the state of self-oscillation is organized, 
and the large-scale oscillation is stabilized. 

We derive the above dynamic self-consistent solution 
on the largest scale. Let g{r) be a function rapidly de- 
creasing and periodic on D, which express the distribu- 
tion of the point vortices. Since the concentrated vortices 
remain inside the 1-to-l resonance islands of tori, the dis- 
tribution function is supposed to be advected along the 
stream line, without changing its form, g, in the first ap- 
proximation. Then the vorticity field attributed to the 
point vortices is, 



■{t,r)-n{r) 



-g{r - R-) - g{r + R.) 
+gir - R+) + g{r + R+) 



(8) 



where ±R± are the coordinates of the four resonance 
islands; R- = {-X,Y) and R+ = {-X + + 
^). In the first approximation, {X{t), Y{t)) moves along 
the base stream line; Y(t) is constant with respect to 
time, Yq, and accordingly X{t) changes in a constant 
speed, X{t) = + 1 sm{2TrYo / Ly) . Then the Fourier 
components of the vorticity field are, 

d;P-^-(i) + gg{-2cos{q^X - QyY) 

+ 2{-l) — i-l)^ COs{q^X + qyY)}. (9) 

In particular, ZP-^-{t) = —Re w^'j^^g) ~ 4g(i q) cosX(t) = 
4^(1,0) sin(t sin(27rFo/I/j,)). Thus the self consistent solu- 
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tion, w"*"^ = cjP '^- , is yielded, on the largest scale, as, 

& = 2g(i,o), — =sin— (10) 



In conclusion, we have discovered the non-stationary 
flow in the 2D Euler equations. The motion of positive 
and negative coherent vortices leads to the large-scale 
oscillation. In the vicinity of the bifurcation point, the 
amplitude of oscillation of the order parameter Z{t) in- 
creases proportionally to (e — Cc)^/^, suggesting Hopf bi- 
furcation. This implies a structures of low-dimensional 
dissipative dynamical system in the high-dimensional 
conservative dynamics of the 2D Euler equations. By 
comparing the dynamics of these coherent vortices in the 
original 2D Euler equations with the one-body dynamics 
of the test point vortex under the auxiliary oscillatory 
fields, we have explained self-consistcntly that this oscil- 
lation result from the collective organization of the state 
of self-oscillation. 

We have focused our consideration only on the largest 
scale. In turbulence, in general, the large-scale and small- 
scale dynamics are not simply separated. Indeed, as 
shown in FIG. 2, neither the oscillatory nor the stationary 
state is completely periodic nor constant, respectively, 
but is accompanied by small fluctuations; accordingly, 
the bifurcation is not completely clear (FIG. 5). These 
small fluctuations are likely to be attributed to the os- 
cillations within the single resonance islands [18] as well 
as the motions of higher-to-higher resonances. Taking 
larger wavenumber or larger frequency modes of field into 
account could reproduce the oscillation more precisely. 
Nevertheless, it is remarkable that still the mechanism of 
the present large-scale oscillation is essentially explained 
only with the large-scale description. 

The collective oscillation discussed in this paper has 
the same mechanism as that in many-body Hamiltonian 
systems [16]. That is, both these two systems have long- 
range interactions [19], which indeed organize the state 
of self-oscillation. Thus this collective oscillation may be 
a universal phenomenon in long-range interacting con- 
servative systems, whether the system is a many-body or 
continuous system. The benefit of the 2D fluid is that it 
is experimentally confirmed more easily than the Hamil- 
tonian systems. Our initial condition will be reaUzed 
as follows; we prepare the system with the aspect ratio 
unity, in which case the field (2, 3, 4) is stable, and then 
we suddenly change the aspect ratio into F. 

Finally, we also observe a similar motion of coherent 
vortices in the 2D stochastic Navier-Stokes (SNS) equa- 
tions with weak Rayleigh friction and viscosity. Starting 
from arbitrary initial conditions, after a transient, the 
system self-organizes the dynamic pattern. Due to the 



stochasticity, this pattern does not last so long, but dis- 
appear in a rather short time. However, it re-appears 
soon. The annihilation and creation of these pairwise 
vortices repeats continually. 
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